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We present a measurement of the unitarity triangle angle 03 using a Dalitz plot analysis of the 
K%-k'^-k~ decay of the neutral D meson from the B -^ jjy*) j<[\*)± process. The method employs the 
interference between Z)° and D" to extract the angle 03, strong phase 5 and the ratio r of suppressed 
and allowed amplitudes. We apply this method to a 357 fb~^ data sample collected by the Belle 
experiment. The analysis uses three modes: B^ — > DK'^ , B^ — > D* K^ with D* -^ Dtt", and B"*" -^ 
DK*^ with K*'^ — > K%-K^ , as well as the corresponding charge-conjugate modes. From a combined 
maximum likelihood fit to the three modes, we obtain 03 = 53° t}^o (stat) ± 3°(syst) ± 9° (model). 
The corresponding two standard deviation interval is 8° < 03 < 111°. 

PACS numbers: 12.15.Hh, 13.25.Hw, 14.40.Nd 

I. INTRODUCTION 

Determinations of the Cabibbo-Kobayashi-Maskawa (CKM) [1] matrix elements provide important checks on the 
consistency of the standard model and ways to search for new physics. The possibihty of observing direct CP violation 
in i? — > DK decays was first discussed by I. Bigi, A. Carter and A. Sanda [2]. Since then, various methods using CP 
violation in B — > DK decays have been proposed [3-6] to measure the unitarity triangle angle 03. These methods 
are based on two key observations: neutral Z?° and D^ mesons can decay to a common final state, and the decay 
B^ —^ DK~^ can produce neutral D mesons of both fiavors via b — > cus and b —^ ucs transitions, with a relative 
phase 0+ between the two interfering amplitudes that is the sum, 5 + 03, of strong and weak interaction phases. For 
the decay B^ -^ DK~, the relative phase is 6'_ = (5 — 03, so both phases can be extracted from measurements of 
such charge conjugate B decay modes. However, the use of branching fractions alone requires additional information 
to obtain 03. This is provided either by determining the branching fractions of decays to flavour eigenstates (GLW 
method [3]) or by using different neutral D final states (ADS method [6]). 

A Dalitz plot analysis of a three-body final state of the D meson allows one to obtain all the information required 
for determination of 03 in a single decay mode. The use of a Dalitz plot analysis for the extraction of 03 was first 
discussed by D. Atwood. I. Dunictz and A. Soni. in the context of the ADS method [6]. This technique uses the 
interference of Cabibbo-favored D'^ -^ K^ti^tt^ and doubly Cabibbo-suppressed Z?" — > A'^tt+tt" decays. However, 
the small rate for the doubly Cabibbo-suppressed decay limits the sensitivity of this technique. 

Three body final states such as AT^tt+tt" [7, 8] have been suggested as promising modes for the extraction of 03. 
In the Wolfcnstcin parameterization of the CKM matrix elements, the weak parts of the amplitudes that contribute 
to the decay B+ -> DK+ are given by V;;i/„, -- AX^ (for the D°K+ final state) and V;*f,V^^ - A\^{p + irj) (for 
D'^K'^). The two amplitudes interfere as the D'^ and D° mesons decay into the same final state K'^ti^tt^; we denote 
the admixed state as D+. Assuming no CP asymmetry in neutral D decays, the amplitude of the D+ decay as a 
function of Dalitz plot variables m^ = rni,o + and m?_ = tti^o - is 

M+ ^ f{ml,m'i) + re''t"'+'^f{mi,ml), (1) 

where f{m?^,rn^) is the amplitude of the D^ -^ A'^tt+tt^ decay, and r is the ratio of the magnitudes of the two 
interfering amplitudes. The value of r is given by the ratio of the CKM matrix elements I^J^V'^sl/lV'^hKisI ^ ^.38 and 
the color suppression factor, and is estimated to be in the range 0.1-0.2 [9]. 
Similarly, the amplitude of the £*_ decay from B~ -^ DK'^ process is 

M_ = f{ml,ml) + re-''l'^+'^f{ml,m:i). (2) 



The 73° — > KgTr~^Tr~ decay amplitude / can be determined from a large sample of flavor-tagged D" -^ Jf^Tr+Tr" 
decays produced in continuum e~^e~ annihilation. Once / is known, a simultaneous fit of B^ and B~ data allows the 
contributions of r, 03 and S to be separated. The method has a two-fold ambiguity: ((^3, S) and (03 + 180°, S + 180°) 
solutions cannot be separated. We always choose the solution with < 03 < 180°. References [7] and [10] give a more 
detailed description of the technique. 

The method described above can be applied to other modes as well as B^ -^ DK^ decay and its charge-conjugate 
mode (charge conjugate states are implied throughout the paper). Excited states of neutral D and K mesons can 
also be used, although the values of 5 and r can differ for these decays. Previously the Belle [10, 11] and BaBar [12] 
collaborations performed analyses using this technique for _B+ — + DK^ and i?+ -^ D*K^ decays. Both Belle [14] 
and BaBar [13] have also performed analyses of the B^ -^ DK*'^ mode, but the Belle result was not combined with 
that from the B~^ — > D^*'^ K^ modes. The latest Belle analyses [11, 14] were based on a 253 fb^^ data sample. In the 
current paper, we report a measurement of 03 with the combination of B'^ — > DK^, i?+ — *■ D*K^ and B~^ — *■ DK*~^ 
modes based on a 357 fb~^ data sample. This analysis supersedes previous Belle results on 03 using Dalitz plot 
analysis of B+ -^ D^*)Ki*)+ decays. 

II. EVENT SELECTION 

We use a 357 fb^^ data sample, corresponding to 386 x 10^ BB pairs, collected by the Belle detector. The decay 
chains B+ -^ DK+, B+ -^ D*K+ with D* -^ Dn" and B+ -^ DK*+ with K*+ -^ A'g7r+ are selected for the analysis; 
the decays B+ ^ Dn+, B+ -^ D*tt+ with D* — > Dtt° and i?° -^ Z3*+7r~ with £)*+ -^ D7r+ serve as control samples. 
The neutral D meson is reconstructed in the KgTr~^Tr~ final state in all cases. We also select decays of D*^ -^ D^ii^ 
produced via the e^e^ -^ cc continuum process as a high-statistics sample to determine the D^ — > A'^tt+tt^ decay 
amplitude. 

The Belle detector is described in detail elsewhere [15, 16]. It is a large-solid-angle magnetic spectrometer consisting 
of a silicon vertex detector (SVD), a 50-layer central drift chamber (CDC) for charged particle tracking and specific 
ionization measurement (dE/dx), an array of aerogel threshold Cerenkov counters (ACC), time-of- flight scintillation 
counters (TOF), and an array of CsI(Tl) crystals for electromagnetic calorimetry (ECL) located inside a supercon- 
ducting solenoid coil that provides a 1.5 T magnetic field. An iron flux return located outside the coil is instrumented 
to detect K^ mesons and identify muons (KLM). 

Charged tracks are required to satisfy criteria based on the quality of the track fit and the distance from the 
interaction point in both longitudinal and transverse planes with respect to the beam axis. To reduce the low 
momentum combinatorial background we require each track to have a transverse momentum greater than 100 MeV/c. 
Separation of kaons and pions is accomplished by combining the responses of the ACC and the TOF with the dE/dx 
measurement from the CDC to form a likelihood C{h) where /i is a pion or a kaon. Charged particles arc identified as 
pions or kaons using the likelihood ratio 7?.pid(/i) = C{h)/{C{K) + C{ti)). For charged kaon identification, we require 
Tl-'PiTiiK) > 0.7. This requirement selects kaons with an efficiency of 80% and pions with an efficiency of 5%. 

Photon candidates are required to have ECL energy greater than 30 MeV. Neutral pion candidates are formed from 
pairs of photons with invariant masses in the range 120 to 150 MeV/c^, i.e. less than two standard deviations from 
the tt" mass. 

Neutral kaons are reconstructed from pairs of oppositely charged tracks without any pion PID requirement. We 
require the reconstructed vertex distance from the interaction point in the plane transverse to the beam axis to be 
more than 1 mm and the invariant mass M^^t^ to satisfy \MTrTT ~ Mj^a\ < 10 MeV/c^, i.e. less than four standard 
deviations from the nominal Kg mass. 

A. Selection of £)*" —y D°tt- 

To determine the D^ -^ X^tt+tt" decay amplitude we use D*^ mesons produced via the e+e" — > cc continuum 
process. The flavor of the neutral D meson is tagged by the charge of the slow pion (which we denote as tt^) in the 
decay!?*- -^D^n-. 

To select neutral D candidates we require the invariant mass of the Kgir^n^ system to be within 9 MeV/c^ of the 
D° mass, M^io. To select events originating from a D*" decay we impose a requirement on the difference AM = 
M^o^+^-^^ — M^o^+^- of the invariant masses of the D*~ and the neutral D candidates: 144.6 MeV/c^ < AM < 

146.4 MeV/c^ The resolutions of the selection variables are ct(AM) = 0.38 MeV/c^ and a(M^o^+^-) = 5.4 MeV/c^ 
The suppression of the combinatorial background from BB events is achieved by requiring the D*^ momentum in 
the center-of-mass (CM) frame to be greater than 2.7 GeV/c. 



The number of events that pass all selection criteria is 271621. To obtain the number of background events 
in our sample we fit the AM distribution. The background is parameterized with the function b{AM) ~ 
(l/AM)^(AA//m7r)2 — 1; the function describing the signal is a combination of two Gaussian peaks with the same 
mean value. The fit finds (261.9 ± 1.1) x lO'^ signal events and 8698 ± 77 background events corresponding to a 
background fraction of 3.2%. 

B. Selection of B+ -» DK+ 

The selection of B candidates is based on the CM energy difference AE = ^Ei — i^boam and the beam-constrained 
B meson mass Mbc = ■\/^bcam ^ (^PiY- where ii^boam is the CM beam energy, and Ei and Pi are the CM energies and 
momenta of the B candidate decay products. We select events with Afbc > 5.2 GeV/c^ and \AE\ < 0.2 GeV for the 
analysis. In addition, we impose a requirement on the invariant mass of the neutral D candidate: |M^o^+^- — A//)o| < 
11 MeV/c2. 

To suppress background from e+e^ ^^ QQ {q = u,d,s,c) continuum events, we require | cos^thrl < 0.8, where ^thr 
is the angle between the thrust axis of the B candidate daughters and that of the rest of the event. For additional 
background rejection, we use a Fisher discriminant composed of 11 parameters [17]: the production angle of the B 
candidate, the angle of the B thrust axis relative to the beam axis and nine parameters representing the momentum 
flow in the event relative to the B thrust axis in the CM frame. We apply a requirement on the Fisher discriminant 
that retains 90% of the signal and rejects 40% of the remaining continuum background. 

The AE and Mbc distributions for _B+ -^ DK~^ candidates are shown in Fig. la,b. The peak in the AE distribution 
at AE = 50 MeV is due to B~^ — > Dn^ decays where the pion is misidentified as a kaon. The 5+ -^ DK^ 
selection efficiency (11%) is determined from a Monte Carlo (MC) simulation. The number of events in the signal 
box (|A£;| < 30 MeV, 5.27 GeV/c^ < Afbc < 5.3 GeV/c^) is 470. 

For the selected events a two-dimensional unbinned maximum likelihood fit in the variables Afbc and AE is per- 
formed. The resulting signal and background density functions are used in the Dalitz plot fit to obtain the event-by- 
event signal to background ratio. To parametrize the shape of the AE — Mbc distribution, we use two-dimensional 
Gaussian peaks for the signal contribution and B^ -^ Dtt^ background. The non-peaking background is parametrized 
by the sum of two components: the product of an empirical shape introduced by ARGUS [18] as a function of Afbc 
and a linear function in AE, and the product of a Gaussian distribution in Mbc and a linear function in AE. Only 
the region with AE > —0.1 GeV is used in the fit; the region with AE' < —0.1 GeV includes a large fraction of events 
from B -^ D*K decay with a lost pion that do not contribute to the signal box. The number of events in the signal 
peak obtained from the fit is 331 ± 23; the event purity in the signal box is 67%. 

C. Selection of B+ -^ D* K+ 

For the selection of i?+ -^ D*K~^ events, in addition to the requirements described above, we require that the mass 
difference AM = Mj(0^+^-^o-MiiO^+^- of neutral D* and D candidates satisfies 140 MeV/c^ < AAf < 145 MeV/c^ 
Figures lc,d show the AE' and Mbc distributions for i?+ -^ D*K^ candidates. The selection efficiency is 6.2%. The 
number of events in the signal box is 111. The parametrization of background and signal shapes is similar to that in 
the 5+ — > DK~^ case. The number of events in the signal peak obtained from the fit is 81 ± 11; the event purity in 
the signal box is 77%. 

D. Selection of B+ -^ DK*+ 

For B^ — > DK*^ decay candidates, in addition to the requirements placed on the B^ -^ DK^ mode, A'*+ 
candidate selection is performed. To select K*^ meson candidates, we require the Kgir^ invariant mass to be within 
50 MeV/c^ of the K*~^ nominal mass. For continuum background suppression, in addition to the |cos0thr| < 0.8 
condition, we require that the K*~^ helicity angle satisfy [cos^hoij > 0.4. The A'*+ helicity angle is defined as the 
angle between the axis of the K*~^ decay products, and the momentum of the B meson, in the K*'^ rest frame. We 
also apply a requirement on the Fisher discriminant that retains 95% of the signal and rejects 30% of the remaining 
continuum background. Figures le,f show the AE and Mbc distributions for B~^ -^ DK*^ candidates. The selection 
efficiency is 4.1%. The number of events in the signal box is 78. The parametrization of background and signal shapes 
is similar to that in the i?+ -^ DK~^ case, but without the contribution of K/n misidentification background. The 
number of events in the signal peak obtained from the fit is 54 ± 8; the event purity in the signal box is 65%. 
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FIG. 1: AE and Mbc distributions for the (a,b) B+ -^ DK+, (c,d) B+ -^ D*K+, and (e,f) B+ -^ DK*+ event samples. 
Points with error bars are the data, and the histogram is the result of a MC simulation according to the fit result. 



III. DETERMINATION OF THE D° -* K^tt+tt- DECAY AMPLITUDE 



s^ 



The amplitude / for the D^ -^ Kgn^ir decay is described by a coherent sum of two-body decay amphtudes and 
one non-resonant decay amphtude, 

N 

f{ml,ml) = ^aje'"^ Aj{ml,ml) + be'^ , (3) 

where N is the total number of resonances, Aj{m'^,m?_) is the matrix element, a^ and aj are the amplitude and 
phase of the matrix element, respectively, of the j-th resonance, and b and /3 are the amplitude and phase of the 
non-resonant component. The total phase and amplitude are arbitrary. To be consistent with other analyses [12, 19] 
we have chosen the D'^ -^ ^sP mode to have unit amplitude and zero relative phase. The description of the matrix 
elements follows Rcf. [20]. 

We use a set of 18 two-body amplitudes. These include five Cabibbo-allowed amplitudes: ii'*(892)+7r^, 
ii:*(1410)+7r-, K^{U30)+7r-, K^{U30)+Tr- and i^*(1680)+7r-; their doubly Cabibbo-suppressed partners; and eight 
amplitudes with K^ and a tttt resonance: K°p, K'^uj, is:^/o(980), Xg/2(1270), i4:^/o(1370), A'0p(1450), K'^ai and 
Kg<J2- The differences from our previous publication [10] are, i) the addition of Kgp{1450), iir*(1410)+7r~ and the 
corresponding doubly Cabibbo-suppressed mode, ii) the use of the Gounaris-Sakurai [21] amplitude description for 
the Kgp and Kgp{l450) contributions, and iii) the mass and width for the /o(1370) state are now taken from Ref. [22] 
{M = 1434 MeV/c2, T = 173 MeV/c^j. 

We use an unbinned maximum likelihood technique to fit the Dalitz plot distribution to the model described by 
Eq. 3. We minimize the negative logarithm of the likelihood function in the form 



-21ogL 



n „ 

4=1 n 



(4) 



where i runs over all selected event candidates, and ?Ti^ j, m}_ ^ are measured Dalitz plot variables. The integral in 
the second term accounts for the overall normalization of the probability density. 
The Dalitz plot density is represented by 

p(TO^,m^) — e(TO^,m^)|/(m^,TO^)|^ + B{m'^,m'^_), (5) 

where f(m']^,m?J) is the decay amplitude described by Eq. 3, e{rn?^,rn?_) is the efficiency, and B(m?^,m'^_) is the 
background density. To take into account the finite momentum resolution of the detector, the Dalitz plot density 
described by Eq. 5 is convolved with a resolution function. The free parameters of the minimization are the amplitudes 
aj and phases aj of the resonances (except for the K^p component, for which the parameters are fixed), the amplitude 
b and phase (3 of the non-resonant component and the masses and widths of the ui and CT2 scalars. 

The procedures for determining the background density, the efficiency, and the resolution of the squared invariant 
mass, are the same as in the previous analyses [10, 11]. The D^ -^ Kgn^n^ Dalitz plot distribution, as well as its 
projections with the fit results superimposed, are shown in Fig. 2. The fit results are given in Table I. The fit fractions 
quoted in Table I for specific modes are defined as the integral of the absolute value squared of the individual mode 
divided by the integral of the absolute value squared of the total amplitude. Due to the interference effects these fit 
fractions may not sum to unity. 

The parameters obtained for the cri resonance (A/ct^ = 519±6 MeV/c^, Fo-i = 454±12 MeV/c^) are similar to those 
observed by other experiments [19, 23]. The second scalar term (72 is introduced to account for structure observed at 
ml^ '- 1.1 GeV^/c^: the fit finds a smaU but significant contribution with M^^ = 1050 ± 8 MeV/c^, T^^ = 101 ± 7 
MeV/c^ (the errors are statistical only). The large peak in the m^ distribution corresponds to the dominant D° -^ 
iir*(892)+7r~ mode. The minimum in the m?. distribution at 0.8 GeV^/c^ is due to destructive interference with the 
doubly Cabibbo suppressed Z?° — > i^*(892)~7r+ amplitude. In the m^^ distribution, the D'^ — > Kgp contribution is 
visible around 0.5 GeV^/c^ with a steep edge on the upper side due to interference with D'^ — > KgUj. The minimum 
around 0.9 GeV^/c^ is due to the decay D^ — > Kgfo{98Q) interfering destructively with other modes. 

We perform a x^ test to check the quality of the fit, dividing the Dalitz plot into square regions 0.05 x 0.05 GeV^/c'*. 
The test finds a reduced chi-square 'X^ jndf = 2.72 for 1081 degrees of freedom {ndf), which is large. Examining Fig. 2, 
we find that the main features of the Dalitz plot are well-reproduced, with some significant but numerically small 
discrepancies at peaks and dips of the distribution. In our final results we include a conservative contribution to the 
systematic error due to uncertainties in the D'^ decay model, discussed in Section IV G. 
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FIG. 2: (a) 7n_, (b) m+, (c) m^^ and (d) Dalitz plot distribution for D* — > D'V^ , D'^ ^ I 
continuum process. The points with error bars show the data; the smooth curve is the fit resuh. 



■K decays from the e^e -^ cc 



IV. DALITZ PLOT ANALYSIS OF B+ -^ £)(*)a:(*'+ DECAYS 



In our previous analyses, the two Dalitz distributions corresponding to the decays of B^ and B^ were fitted 
simultaneously to give the parameters r, (jj^ and d. Confidence intervals were then calculated using a frequentist 
technique, relying on toy MC simulation. In this approach, there was a bias in the fitted value of the (positive 
definite) parameter r, and the errors on 03 and 6 were also r-dependent. 

In the present analysis, we use a method similar to that of BaBar [12]: fitting the Dalitz distributions of the B'^ 
and B~ samples separately, using Cartesian parameters x± — r± cos{±(f>T, + 6) and y± = r± sin(±03 + 6), where the 
indices "+" and "— " correspond to B+ and B~ decays, respectively. Note that in this approach the amplitude ratios 
(r_|_ and r_) are not constrained to be equal for the B~^ and B~ samples. Confidence intervals in r, 03 and S are 
then obtained from the {x±, y±) using a frequentist technique. The advantage of this approach is low bias and simple 
distributions of the fitted parameters, at the price of fitting in a space with higher dimensionality {x+,y^,X-,y-) 
than that of the physical parameters (r, 4>3,S); see Section IV E. 

The fit to a single Dalitz distribution with free parameters x and y is performed by minimizing the negative unbinncd 



TABLE I: Fit results for D" — > K'^tv^tt decay. Errors are statistical only. 



Intermediate state 


Amplitude 


Phase (°) 


Fit fraction 


T^gai 


1.43 ±0.07 


212 ±3 


9.8% 


K'sP° 


1.0 (fixed) 


(fixed) 


21.6% 


K°suj 


0.0314 ± 0.0008 


110.8 ±1.6 


0.4% 


K°sM980) 


0.365 ± 0.006 


201.9 ±1.9 


4.9% 


K°a2 


0.23 ± 0.02 


237 ±11 


0.6% 


K°s f2{1270) 


1.32 ±0.04 


348 ±2 


1.5% 


K°sMl370) 


1.44 ±0.10 


82 ±6 


1.1% 


K°sP°{U50) 


0.66 ± 0.07 


9±8 


0.4% 


K*{892) + n'- 


1.644 ±0.010 


132.1 ±0.5 


61.2% 


K''{892)-n+ 


0.144 ±0.004 


320.3 ±1.5 


0.55% 


Ji'*(1410) + 7r- 


0.61 ±0.06 


113 ±4 


0.05% 


ii'*(1410)-7r+ 


0.45 ± 0.04 


254 ±5 


0.14% 


is:o*(1430) + 7r- 


2.15 ±0.04 


353.6 ±1.2 


7.4% 


ii:o(1430)"7r+ 


0.47 ± 0.04 


88 ±4 


0.43% 


is:|(1430) + 7r- 


0.88 ± 0.03 


318.7 ±1.9 


2.2% 


A'2*(1430)-7r+ 


0.25 ± 0.02 


265 ±6 


0.09% 


J\*(1680) + 7r- 


1.39 ±0.27 


103 ±12 


0.36% 


J-i:*(1680)"7r+ 


1.2 ±0.2 


118±11 


0.11% 


non-resonant 


3.0 ±0.3 


164 ±5 


9.7% 



likelihood function 



-21ogL = -2 



4=1 



logp(m+^„ ml,, AE^, Afbc,*) - log / p{m+,mi,AE, Mbc)) dmt^ dmt dAE dM^c 



D 



(6) 



with the Dalitz distribution density p represented as 

p{ml, ml,AE, A/bc) = e{ml,ml)\f{ml, ml) + {x + iy)f{ml,ml)\^F,ig{AE, M^,) + F\,,^{ml,ml, AE, A/bc), (7) 

where i^sig is the signal distribution as a function of AE and Afbc (represented by the product of two Gaussian shapes), 
i'bck is the distribution of the background, and e(m']^^m?_) is the efficiency distribution over the phase space. As in 
the study of the sample from continuum D*~ — > D^ttJ decays, the finite momentum resolution is taken into account 
by convolving the function (7) with a Gaussian resolution function. The efficiency and the momentum resolution were 
extracted from the signal MC sample, where the neutral D meson decays according to phase space. The determination 
of the background contribution and efficiency profile is described below. 



A. Backgrounds 



To take backgrounds into account in the analysis, their Dalitz plot distributions, AE — Mbc distributions (which 
in general may depend on Dalitz plot region) and relative fractions have to be known. The backgrounds are divided 
into three categories: 



Continuum background 



• B^ -^ I?'*V+ background (with tt /K misidcntification) . This background is relevant only for the 5+ 
and B+ -^ D*K+ modes. 



DK^ 



• Other hh backgrounds. 

Continuum background (from the process e+e~ — > gg, where q = u, d, s, c) gives the largest contribution. It 
includes both pure combinatorial background, and continuum Z?° mesons combined with a random kaon. This type of 
background is studied in data with cos ^thr and Fisher discriminant requirements applied to select continuum events. 



TABLE II: Results of fits to test samples in parameters {x,y). Errors are statistical only. 



Mode 


x_ 


y- 


x+ 


y+ 


B+ ^ Dt:+ 


-0.030 ±0.015 


-0.014 ±0.016 


-0.003 ±0.014 


-0.021 ±0.018 


B+ -^ D'n+ 


-0.007 ± 0.035 


-0.078 ± 0.039 


0.004 ± 0.032 


0.022 ± 0.039 


B ^D*--K+ 


0.016 ± 0.030 


0.014 ± 0.030 


0.026 ± 0.026 


0.009 ± 0.036 



To check that the Dalitz plot shape of the continuum background selected using these requirements corresponds to 
that of the signal region, we use a MC sample that includes e"'"e~ —>■ qq {q = u,d,s,c) decays. The Dalitz plot 
distribution of the continuum background is parametrized by a third-order polynomial in the variables rn^ and m'^ 
(which represents the combinatorial component) and a sum of D'^ and D'^ shapes for real neutral D mesons combined 
with random kaons. 

The B^ -^ D^*>tt'^ process with a pion misidentified as a kaon is suppressed by the requirements on the A'/tt 
identification variable T^pid and the CM energy difference. The fraction of this background is obtained by fitting 
the AiJ — Mbc distribution; the corresponding Dalitz plot distribution is that of a Z?" without the opposite flavor 
admixture. 

Other BB backgrounds, of which the dominant fraction comes from the decay of £)(*)" from one B meson, with 
some particles taken from the other B decay, are investigated with e+e^ -^ T(45') — > BB MC events. The Dalitz 
plot distribution of the 66 background is parameterized by a second-order polynomial (for the i?+ -^ DK^ and 



_B+ — > DK*^ modes) or by a linear function (for the B^ 
shape (U" for B+ data and L»° for B' data). 



D*K^ mode) of m^ and to_, plus a "correct-flavor" D^ 



B. Efficiency 

Knowledge of the absolute value of the reconstruction efficiency is not essential for our analysis. However, the 
relative variations of the efficiency over the D^ decay phase space can affect the fit result. The shapes of the efficiency 
across the phase space are studied for each mode using signal MC samples with a constant i?° — > ifg7r+7r~ decay 
amplitude. The efficiency profiles are fitted with third-order polynomial functions of mJi^ and m?_ symmetric under 
exchange of 7r+ and 7r~. The efficiency is nearly fiat over the phase space, falling by 10-20% at its edges (relative to 
the efficiency at the center). 



C. Control sample fits 



To test the consistency of the fitting procedure, the same procedure was applied to the B^ 



b^f''n+ and B 



D* 



'7r+ control samples as to the i?+ -^ Dj^ K'^*^^ signal. For decays to which only one D fiavor can contribute, the 

D\ 7r+ a small amplitude 



fit should return values of the amplitude ratio r consistent with zero. In the case of B^ 
ratio r ~ 0.01 is expected (due to the small ratio of the weak coefficients |^(,V^dl/|V^{,V^„*dl ^ 0.02 and the additional 
color suppression factor as in the case of i?+ -^ D^K^). Deviations from these values can appear if the Dalitz plot 
distribution is not well described by the fit model. 

For the control sample fits, we treat i?+ and B~ data separately, to check for the absence of CP violation. The 
free parameters of the Dalitz plot fit are x± and y± . 

The fit results for the three test samples are presented in Table II; contour plots showing integer multiples of the 
standard deviation in the x± and y± variables for the three test samples are shown in Fig. 3. The results are consistent 
with r ^ 0.01 for the B+ -^ Dtt^ and i3+ -^ D*it^ modes and with zero for the B -^ D*^7r+ mode. 



D. Signal fit results 



The Dalitz distributions of B+ -^ D+K^, _B+ -^ DZ^K^ , and i?+ —f D+K*^ modes for the events in the signal 
region used in the fits are shown in Fig. 4. The results of the separate B^ and B~ data fits are shown in Fig. 5. 
The plots show the constraints on the parameters x and y for the B^ -^ D+K^ , B^ -^ D*^K^ and B'^ — s- D^K*^ 
samples. The values of the fit parameters x± and y± are listed in Table III. 

The fit to the 5+ — » DK^ sample yields r+ and r_ values close to each other (r+ = 0.172 and r_ = 0.160) 
with significantly different values of the total phase 6+ and 6*-, which is an indication of CP violation. For the 
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FIG. 3: Results of test sample fits with free parameters x — rcos6 and y = r sin6' for (a) B'^ —^ Dn'^ , (b) _B^ -^ D*tx^ and 
(c) B — > D*~-k'^ samples, separately for B~ and B"*" data. Contours indicate integer multiples of the standard deviation. 



TABLE III; Results of the signal fits in parameters (a, y). Errors are statistical only. 



Mode 


a;„ 


y- 


x+ 


y+ 


B+ -^ DK+ 
B+ -^ D*K+ 
B+ -^ DK*+ 


U.U/0_Q QgQ 

-0.1281°:}!^ 
-0.784l°:^« 


0.1701?-- 
-0.3391°:}^^ 
-0.2811°:^^° 


-0.135+-.«« 

0.032«;}?» 

-0.105«:}^? 


-0.085«:«^» 

0.008«;1^^ 

-O.OOl^/J^ 



5+ -^ D*K^ and B+ — > DK*^ samples, the values of r_|. and r_ differ (with large errors); in both cases, r+ is within 
a standard deviation of zero, so the corresponding phase ^+ is poorly determined and the significance of CP violation 
is low. 



E. Evaluation of statistical error 



We use a frcqucntist technique to evaluate the statistical significance of the measurements. This method requires 
knowledge of the probability density function (PDF) of the reconstructed parameters x and y as a function of the true 
parameters x and y. To obtain this PDF, we employ a "toy" MC technique that uses a simplified MC simulation of the 
experiment which incorporates the same efficiencies, resolution and backgrounds as used in the fit to the experimental 
data. 

For each mode, 4000 toy MC experiments arc made with generated {x^y) equal to the experimental fit results. 
The PDFs obtained are approximately Gaussian, but we observe the width of the distribution to vary with x and y. 
There is also an offset in the parameters x, y introduced by the fit procedure. The correlations between x and y are 
negligible. For the _B+ -^ DK^ and i?+ -^ D*K^ modes, the deviations from the Gaussian shape are small, and 
taking into account the tails of the distributions has only a minor effect on the fit result. However, the result of the 
measurement in the B^ -^ DK*^ mode has a relatively large value of r (r_ — 0.83 ± 0.28). For such a large value, 
the toy MC study shows the existence of significant tails in the observed x and y distributions. To parametrize the 
shapes of the x and y distributions for the B^ —^ DK*~^ sample we perform additional toy MC simulation for different 
generated x and y, and fit the distributions with the sum of two Gaussians ("narrow" and "wide"). To obtain the 
PDF shape for any given set of {x,y), linear interpolation in x and y (for the bias of the wide Gaussian relative to 
the narrow one) and in x'^ and y^ (for the width and fraction of the wide Gaussian) are used. 

To obtain the (x, y) dependence of a^ and <7y and the offset in x and y, we generate toy MC experiments at different 
points in the {x, y) plane, corresponding to 4 values of 9 for r = 0.15 and 8 values of for r = 0.3 and r = 0.6. For 
the B^ -^ DK*^ mode, an additional set of 8 points with r = 0.9 is added. 200 toy MC experiments are generated 
at each point. The standard deviation and offset are then parametrized with polynomial functions in x and y. 

The offset in x and y is large enough to affect the result of the measurement (up to 0.02 for the B^ -^ DK^ mode); 
it is larger in the i?+ — > D*K^ and i?+ -^ DK*^ modes, where the data samples are smaller. To test the effect of 
sample size, toy MC experiments were generated with event samples ten times larger than those of the B^ -^ DK^ 
data: the offset for this sample is less than 5 x 10~^. This is consistent with the maximum likelihood fit being unbiased 
in the limit of large sample sizes. For each mode, the offset is included in the PDF for the reconstructed parameters, 
so the confidence intervals (and central values) for (jy^, r, S are unbiased. 
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FIG. 4: Dalitz plots for the neutral D meson from B and B^ decays, respectively, for the (a,b) B"*" -^ DK^ , (c,d) B^ 
D'K'^ , and (e,f) B'^ — > DK*^ modes. Note that the axes are flipped in the case of B~ with respect to B^ . 
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FIG. 5: Results of signal fits with free parameters x — rcosd and y — rsinO for (a) B+ — > DK^ , (b) B+ -^ D* K^ and (c) 
B"^ -^ DK*^ samples, separately for B~ and B^ data. Contours indicate integer multiples of the standard deviation. 



TABLE IV: CP fit results. The error intervals are statistical only. 



Parameter 


B+ ^ DK+ mode 


B+ -^ D''K+ mode 


B+ -^ DK*+ mode 


la interval 


2a interval 


Icr interval 


2a interval 


la interval 


2cr interval 


<t>3 

r 
5 


19 9° 

0.165l»:°^^ 

147.5° +--: 


20.8° < 03 < 108.4° 

0.048 < r < 0.281 

100.3° <5 < 187.7° 


86.1° l^Jfo 
0.207l»:-3l^ 

QO^ 7° +36.8° 
OO'i.l -93.6° 


r < 0.510 


10.8° ifj: 

0.5491°;?^^ 
210.2° +f, f. 


0.242 < r < 1.120 



Different implementations of the frequentist technique can be used to obtain the values of the physical parameters. 
The two most widely used arc central (so-called Neyman) intervals, and the unified approach of Fcldman and Cousins 
[24]. The latter takes unphysical regions of the parameter space into account. In this analysis, fitting with four 
parameters {x±,y±) in a problem defined by three parameters (r, 03, 5) creates an unphysical region: all {x±,y±) 
values where r+ 7^ r_. For two of the three samples (B+ — > D*K^ and B^ — * DK*^) the values of r+ and ?-_ differ 
significantly, i.e. the result for these modes is well into the unphysical region. Central intervals would overestimate 
the significance of the measurement in this case. Therefore, to obtain the values of the physical parameters r, 03 and 
5 (we will refer to this set of parameters as a vector /i) given the measurement result a;+, j/+, .t_ and y- (or vector 
z) we use the Feldman- Cousins approach. 

The confidence level a is calculated as 



a{tj) 



p{z\^)dz, 



(8) 



•D(p) 



where the integration domain T) is given by likelihood ratio ordering: 



> 



p{z\fii,astiz)) p{zo\fJ.hcstizo))' 



(9) 



Here, p{z\fj.) is the normalized probability density to obtain the measurement result z for a given set of true parameters 
fi. This PDF is obtained by toy MC simulation, /-ibcst(^) stands for the best true parameters for a given measurement 
z, i.e. ^ such that p{z\fi) is maximized for the given z. zq is the result of the fit to experimental data. 

Fig 6 shows the projections of the three-dimensional confidence regions onto (r, 03) and (03, i5) planes. We show the 
20%, 74% and 97% confidence level regions, which correspond to one, two, and three standard deviations for a three- 
dimensional Gaussian distribution. The central values for the parameters r, 03 and S with their one and two standard 
deviation intervals are presented in Table IV. The confidence intervals in Table IV are calculated from projections of 
the three-dimensional confidence regions onto each of the parameters. Toy MC studies of the B^ —^ DK^ mode show 
that these intervals include the true value of 03 slightly more often than the confidence levels imply (overcoverage) , 
and are therefore conservative. 
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FIG. 6: Projections of confidence regions for the (a,b) B^ —^ DK^ , (c,d) B"*" — > D* K^ , and (e,f) B^ — > DK*^ modes onto 
the (r,(f)-i) and ((jia,!?) planes. Contours indicate integer multiples of the standard deviation. 



TABLE V: Systematic errors in x, y variables. 



Component 


B+ -> DK+ mode 


B+ -^ D''K+ mode 


B+ -^ DK''+ mode 


Ax- 


Ay- 


Aa::+ 


Ay+ 


Aa::- 


Ay- 


Ax+ 


Ay+ 


Ax- 


Ay_ 


Ax+ 


Aj/+ 


Background shape 

Efficiency shape 

AE — Mhc Stat, uncertainty 

AE — Mbc difference for qq and bb 


0.006 
0.011 
0.002 
0.004 


0.015 
0.004 
0.005 
0.002 


0.002 
0.017 
0.002 
0.000 


0.006 
0.001 
0.004 
0.005 


0.011 
0.001 
0.005 
0.019 


0.023 
0.001 
0.009 
0.012 


0.002 
0.001 
0.002 
0.002 


0.004 
0.001 
0.003 
0.010 


0.011 
0.024 
0.006 
0.010 


0.026 
0.023 
0.013 
0.027 


0.002 
0.002 
0.003 
0.004 


0.009 
0.009 
0.003 
0.001 


Total 


0.013 


0.016 


0.017 


0.009 


0.023 


0.027 


0.004 


0.011 


0.029 


0.046 


0.006 


0.013 



F. Estimation of systematic errors 
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Experimental systematic errors come from uncertainty in the knowledge of the functions used in the signal Dalitz 
plot fit. These include the Dalitz plot profiles of the backgrounds and the detection efficiency, the momentum resolution 
description, and, since the background fraction varies depending on the AE and A/bc parameters, the parametrizations 
of the AE — Mbc shape of the signal and background. 

The systematic errors related to the Dalitz plot shape of the background are studied by using different background 
descriptions in the signal fit. Instead of the baseline background shape, in which the continuum background is 
determined from the continuum-enriched experimental sample, and the bb background from a MC sample, we perform 
fits with the background shape represented by only the continuum background parametrization, and the background 
shape extracted from the Mbc sideband (5.2 GeV/c^ < Mbc < 5.26 GeV /(?). The maximum deviation of the 
parameters x, y from those obtained with the "standard" background parameterization was taken as a measure of the 
systematic error. 

The distribution of the relative detection efficiency over the Dalitz plot is obtained from the MC simulation. To 
estimate the effect of imperfect simulation of the detector response, we also attempt to obtain the efficiency shape from 
the experimental data. We expect that the detection efficiency depends primarily on the momentum of the D meson. 
We also assume that for high momentum D mesons, the efficiency does not vary significantly over phase space and is 
well determined in the simulation. Therefore, to obtain the relative efficiency distribution for, e.g. the i?"*" — *■ DK~^ 
sample, we use the continuum D*~ -^ D'^t^7 sample with £)" — > KgTT^TT~ (the one used for the determination of 
the D° — > i^^TT+TT" amplitude) and compare Dalitz plots for high momentum D'^ {po > 3.5 GeV/c) and 1)° in the 
momentum range corresponding to the B+ -^ DK~^ decay (2.1 GeV/c < po < 2.4 GeV/c). If the background is 
subtracted, the ratio of the event densities for the decays of Z)° in different momentum ranges should equal the ratio of 
detection efficiencies. Once we assume that the efficiency for high momentum D'^ is known, the efficiency distribution 
of Z)° from B+ -^ DK^ decay can be calculated. The difference of the parameters x, y between fits using MC-based 
and data-based efficiency shapes is taken as the corresponding systematic error. 

The uncertainty due to momentum resolution is estimated by performing a D" model fit with a subsequent fit to 
B decay data without taking the momentum resolution into account. The bias of the {x, y) parameters is negligible, 
therefore, this systematic uncertainty is neglected. 

The uncertainty of the AE — Mbc distribution parameterization (which, in particular, leads to uncertainty in 
the signal-to-background ratio) is taken into account by performing fits using AE — 7\fbc shapes with each of the 
parameters varied by one standard deviation. 

The AS — Mbc shape and the total background fraction are extracted from the signal AS — Mbc fit. Both continuum 
and bb backgrounds are parametrized by a common AE — Afbc density. However, in general, continuum and hh (and 
even each of the sub-components of these backgrounds) may be distributed differently, and the distribution may be 
correlated with the Dalitz plot distribution. To estimate the contribution of this effect, we make additional fits to 
data with individual AE — Afbc distributions for continuum and bb background components extracted from generic 
MC samples. The corresponding bias is included in the background systematic error. 

The results of the study of experimental systematic errors are summarized in Table V for each mode. The un- 
certainties in the parameters (x, y) are then used to obtain the systematic errors for r, 03 and 5 in the combined 
measurement (see Section IV H). In our previous analyses [10, 11] the systematic error was dominated by a term due 
to the bias in the fit of (03, r, 5) to the 5+ — > Dtt"*" control sample (see Ref. [10] Section V.B). In the present analysis, 
biases in the fit of {x±^y±) to data are small, and are taken into account when forming confidence intervals in 03, 
r, and 5 (section IV E). Fitted values for the B^ -^ D^*'>tt'^ control samples are consistent with expectations (sec- 
tion IV C), with statistical precision comparable to the systematic errors listed in Table V. No additional systematic 
error is assigned. 



TABLE VI: Estimation of the D° 



decay model uncertainty. 



Fit model 


(Ar)max 


(A,^3)max (°) 


(A5)^ax (°) 


F, = Fd = 1 




0.01 


3.1 


3.3 


r(g^) — constant 




0.02 


4.7 


9.0 


K'{892)+, p,uj, A'*(892)- 


, /o(980), A:o(1430), non-res. 


0.05 


8.5 


22.9 


No (Tl 




0.01 


2.6 


4.3 


No (72 




0.01 


0.6 


0.7 


CLEO model 




0.02 


5.7 


8.7 



G. Estimation of model uncertainty 
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The model used for the D^ — > K 



'0. 



decay amplitude is one of the main sources of systematic error for our 



analysis. The amplitude is a result of the fit to the experimental Dalitz plot, however since the density of the plot 
is proportional to the absolute value squared of the decay amplitude, the phase of the complex amplitude is not 
directly measured. The phase variations across the Dalitz plot arc therefore a function of model assumptions and 
their uncertainties may affect the D Dalitz plot fit. 

We use a MC simulation to estimate the effects of the model uncertainties. Event samples are generated according 
to the Dalitz distribution described by the amplitude given by Eq. 7 with the resonance parameters extracted from 
our fit to continuum D° — > Kgir^ir^ data, but to fit this distribution different models for /(m^,r7T,?.) are used. We 
scan the phases cjj^ and S in their physical regions and take the maximum deviations of the fit parameters ((Ar)f[iaxj 
(A(/)3)ijiaxj and (Ai5)jjiax) as model uncertainty estimates. These studies use the value r = 0.15, close to the one 
obtained in the dominant 5+ -^ DK^ mode. 

All the fit models are based on Breit-Wigner parameterizations of resonances as in our default model, but with 
differences in the treatement of the broad components. By default, we use Blatt-Weisskopf form factors for the Z?" 
meson (Fd) and intermediate resonances (F^) and a q^ dependence of the resonance width F (see [20]): to estimate 
the effect of the corresponding theoretical uncertainties, we use a fit model without form-factors, and with constant 
resonant widths F^. We also use a model containing only the largest doubly-Cabibbo-suppressed term K* {^%2)~ tt^ , 
the narrow resonances, iCg/o(980), and -K'g(1430)'''7r~, with the remainder of the amplitude approximated by the flat 
non-resonant term. 

Other models used are the model with all the resonances from Tabic I excluding the cri or 1J2 states, and the model 
used by CLEO [19]. The results of the study of model uncertainty are summarized in Table VL The maximum 
deviation of the parameters is taken as the uncertainty due to the Z?" -^ A'g7r"*"7r~ decay model. Each of the 
models listed in Table VI has a reduced y^ significantly poorer than that of the default model. Therefore, our 
model uncertainty does not include an additional contribution due to the poor quality of the model fit, discussed in 
Section IIL 

The mode i?+ — > DK*^ has an additional uncertainty due to the possible presence of a nonresonant B^ -^ DKgiT^ 
component, which can also be treated as a model uncertainty. Since the nonresonant decay is described by the same 
set of diagrams as B'^ -^ DK*^, a similar CP violating effect should be present that in general has values of r and 5 
which differ from those for the resonant mode. Thus, for a (^3 measurement without taking the B^ — > DKgir^ mode 
into account, its contribution can bias the fit parameters. To estimate the corresponding systematic uncertainty, we 
first measure the fraction of nonresonant decays within the i?+ -^ DK*^ signal region. To increase the sample size 



for this study, we include additional D'^ decay modes: K~^Tr 



K'^n+TT- 



K'^TT 7r° and K~^Tr 



Based on the S^ 



yield in the K*~^ mass sidebands, and the observed shape of the invariant mass M^o^+ of signal candidates, we find an 
upper limit of 6.3% on the B^ -^ DKgiT^ fraction. Wc then perform a toy MC simulation with a 6.3% nonresonant 
contribution added to determine the bias of the fit parameters. The fits are performed for various values of the r 
and 5 parameters of the nonresonant component and various values of the relative phase between B — s- DK and 
B^ -^ DKgTT^ amplitudes. The maximum bias of the fit parameters is taken as the corresponding systematic error: 
Ar — 0.084, A(/)3 = 8.3°, AS = 49.3°. The ^3 bias is significantly smaller than that for the strong phase S, since 03 
is obtained from a difference of the total phases for i?+ and B~ decays, and a part of the bias cancels in this case. 



H. Combined 03 measurement 



The three event samples, i?+ -^ DK^, B+ — > D*K^, and S+ -^ DK*~^ are combined in order to improve the 
sensitivity to 03 . The confidence levels for the combination of three modes are obtained using the same technique as 
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TABLE VII: Results of the combination of B^ 



DK+, B^ 



D''K+., and B^ 



DK* modes. 



Parameter 


la statistical interval 


2cr statistical interval 


Systematic error 


Model uncertainty 


03 


rq qO +14.8° 


11.7° <(t>z< 107.7° 


2.5° 


8.7° 


TDK 


n 1 c:n+0-054 

u.ioy_o.050 


0.048 < TDK < 0.271 


0.012 


0.049 


Sdk 


145 7° +19-0° 


100.6° < 5dk < 185.9° 


3.0° 


22.9° 


ro'K 


n ■I7C+0.108 

U-J^ '0-0.099 


< rn'K < 0.407 


0.013 


0.049 


Sd'k 


302.0° +f^% 


- 


6.1° 


22.9° 


TDK' 


0.564l»:?^«5 


0.231 < r-DK' < 1.106 


0.041 


0.084 


Sdk- 


242 6° +^''-^° 

Z1Z.U _23.2o 


186.0° < Sdk" < 300.2° 


2.5° 


49.3° 



for the single mode. In this case, the vector of physical parameters /i — (03, roK-, Sdk, fD*K, Sd*k, tdk*, ^dk*) and 
there are twelve measured parameters: four parameters {x±, y±) for each of the three modes. The probability density 
of the measurement result p{z\n) is the product of the probability densities for the individual modes. 

For each /x, to obtain the confidence level a{ii), the integral in Eq. 8 is performed over a 12-dimensional space. This 
requires extensive computation that makes it impractical to scan the whole range of physical parameters to obtain 
multidimensional confidence regions. However, it is still possible to calculate the confidence intervals for each of the 
individual parameters. 

To calculate the systematic errors for the combined measurement, we vary the measured parameters x and y within 
their systematic errors. Gaussian distributions are used for the variation, and the systematic biases are assumed to 
be uncorrelated. For each of the varied parameter sets, the central values of the physical parameters are calculated. 
The systematic error of the physical parameter is then taken to be equal to the RMS of the resulting distribution. 

The error due to the uncertainty in the D'' -^ /v^tt+tt" amplitude is considered to be equal in all modes. The 
uncertainty due to the possible contribution of the nonresonant 5+ — > DKg7T~^ amplitude to B'^ — > DK*~^ decay 
is also included in the model error, by varying the parameters {x±,y±) of the B'^ —>■ DK*~^ sample according to 
the biases Ar, Acfi^, and A6 obtained above (Section IV G). We find a contribution of only 1.8° to the combined 
03 measurement; this is added in quadrature, yielding a total error of 8.7°. For the r and 6 parameters of the 



S^ 



DK*+ mode the B+ ^ DK9 



contribution dominates the model error. 



Confidence intervals for the combined measurement together with systematic and model errors are shown in Ta- 
ble VII. The confidence intervals are statistical only and are calculated from projections of the seven-dimensional 
confidence regions onto each of the parameters. The statistical confidence level of CP violation is 74%. It is defined as 
the minimum confidence level a{p) for the CP-conserving set of physical parameters ^ {i.e. the set with 03 = 0). The 
two standard deviation interval (including model and systematic uncertainties) is calculated by adding the doubled 
systematic and model errors in quadrature to the two standard deviation statistical errors: we find 8° < 03 < 111°. 

We note that the two standard deviation interval of 03 is more than twice as wide as the one standard deviation 
interval. The non-Gaussian errors in 03 are related to the low significance of r; this effect should be reduced with 
larger data samples in the future. There is also a contribution from the disagreement between the B^ — > DK~^ and 
_B+ — > DK*~^ results, which leads to a second minimum in a(/i) between the one and two standard deviation levels. 



CONCLUSION 



We report the results of a measurement of the unitarity triangle angle 03, using a method based on Dalitz plot 



analysis of Z?" -^ K 



'0. 



decay in the process B^ -^ ]j{*)x^*^~^. The measurement of 03 using this technique was 



performed based on a 357 fb ^ data sample collected by the Belle detector. From the combination of B^ -^ DK^ , 
B+ -^ D*K+ and B+ -^ DK*+ modes, we obtain the value 03 = 53° l^gl ± 3° ± 9°; of the two possible solutions 
we choose the one with < 03 < 180°. The first error is statistical, the second is experimental systematics and the 
third is model uncertainty. The two standard deviation interval (including model and systematic uncertainties) is 
8° < 03 < 111°. The statistical significance of CP violation for the combined measurement is 74%. The method 
allows us to obtain a value of the ratio of the two interfering D decay amplitudes r, which can be used in other 03 
measurements. We obtain r = 0.159taJ^5o ± 0.012 ± 0.049 for the B+ -^ DK+ mode, r = 0.175l^;J[^^ : 
for the B+ -> D*K+ mode and r = 0.564l°;?i6 ± 0.041 ± 0.084 for the B+ -> DK*+ mode. 



17 



Acknowledgments 



We thank the KEKB group for the excellent operation of the accelerator, the KEK cryogenics group for the efficient 
operation of the solenoid, and the KEK computer group and the National Institute of Informatics for valuable 
computing and Super-SINET network support. We acknowledge support from the Ministry of Education, Culture, 
Sports, Science, and Technology of Japan and the Japan Society for the Promotion of Science; the Australian Research 
Council and the Australian Department of Education, Science and Training; the National Science Foundation of China 
and the Knowledge Innovation Program of Chinese Academy of Sciencies under contract No. 10575109 and IHEP-U- 
503; the Department of Science and Technology of India; the BK21 program of the Ministry of Education of Korea, 
and the CHEP SRC program and Basic Research program (grant No. RGl-2005-000-10089-0) of the Korea Science 
and Engineering Foundation; the Polish State Committee for Scientific Research under contract No. 2P03B 01324; 
the Ministry of Science and Technology of the Russian Federation; the Slovenian Research Agency; the Swiss National 
Science Foundation; the National Science Council and the Ministry of Education of Taiwan; and the U.S. Department 
of Energy. 



APPENDIX: MODEL UNCERTAINTY IN (x, y) PARAMETERS 



Averaging the 03 measurements between different methods and different experiments is most conveniently 
done using the {x,y) parameters, which are better-behaved than (r, 03,i5). This requires errors (includ- 
ing model uncertainty) to be expressed in terms of the {x,y) parameters, but our model uncertainty study 



was performed in (r, 4>3,S). To calculate the errors due to the D" 



O^+n 



KU 



decay model uncer- 



tainty in (a;,y) we propagate the (r, 03,(5) errors under the assumption that there is no correlation between 
them. This approach results in substantial off-diagonal terms. The correlation matrix for the vector z = 
{xdk- 1 Vdk- , xdk+ , yDK+ , xd'k- , Vd'K- , xu'K+ 1 yD-K+ , xdk— , Vdk— , xdk*+ , yDK*+ ) is defined as 



and is equal to 



/ 



Cd" ~ 



\ 



459 


10 


10 


241 


104 


-232 


332 


-64 


-56 


-3 


21 


1 


6 





62 


2 


-12 





68 


3 


64 


3 


30 


2 



-104 332 

-232 -64 

264 -68 

-68 436 

-18 55 

7 -22 

2 -4 

20 -57 

-3 9 

22 -66 

20 -59 

10 -30 



(C_Do).j 



-56 21 

-3 1 

-18 7 

55 -22 

515 -108 

-108 282 

-54 -235 

400 -135 



AzjA: 



■3' 



10 
-70 
-64 
-33 



-4 
27 
25 
13 



6 



2 

-4 

-54 

-235 

242 

26 

-1 

7 

5 

2 



62 

2 

20 

-57 

400 

-135 

26 

556 

-13 

76 

69 

33 



-12 68 64 30 \ 

3 3 2 

-3 22 20 10 

9 -66 -59 -30 

10 -70 -64 -33 

-4 27 25 13 

-17 5 2 

-13 76 69 33 

939 36 -328 878 

36 732 -556 -99 

-328 -556 769 -86 

878 -99 -86 904/ 



X 10" 



(A.l) 



The uncertainty due to the non-resonant contribution to the B^ —>■ DK*^ mode produces an additional contribution 



C 



DK-R 



( 766 -371 113 828 \ 
-371 2984 -2605 -1136 

113 -2605 2601 919 

V 828 -1136 919 1150/ 

to the terms in {xdk— ,yDK— iXdk-+ iyDK'+)- 



X 10" 



(A.2) 
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